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Abstract 

We discuss meson-meson scattering at next-to-next-to-leading order in the 
chiral expansion for QCD-like theories with general n degenerate flavours for 
the cases with a complex, real and pseudo-real representation. I.e. with global 
symmetry and breaking pattern SU{n)L x SU{n)R — )• SU{n)v, SU{2n) — )• 
SO{2n) and SU{2n) — )• Sp{2n). We obtain fully analytical expressions for 
all these cases. We discuss the general structure of the amplitude and the 
structure of the possible intermediate channels for all three cases. We derive 
the expressions for the lowest partial wave scattering length in each channel and 
present some representative numerical results. We also show various relations 
between the different cases in the limit of large n. 



PACS: 



1 Introduction 



In an earlier paper [1] we started the phenomenology of QCD-like theories at next- 
to-next-to-leading (NNLO) order in the light mass expansion in their respective low- 
energy effective theories. The motivation for this work is that these theories are 
interesting as variations on QCD and could play some role as models for a nonper- 
turbative Higgs sector. Early work in this context are the technicolor variations of 
[2, 3, 4]. Recent reviews of more modern developments are [5, 6]. Lattice calculations 
have started to explore these type of theories as well, some references are [7]. The 
main interest in these theories is in the massless limit but lattice simulations are nec- 
essarily performed at a finite fermion mass. In [1] we worked out a number of simple 
observables, the mass, decay constant and vacuum-expectation-value to NNLO in 
these theories. Here we work out the amplitude for meson-meson scattering to the 
same order. In lattice calculations the amplitude for meson-meson scattering is not 
directly accessible but the scattering lengths can be derived from the dependence on 
the volume of the lattice [8]. We therefore also provide explicit expressions for the 
scattering lengths. 

The EFT relevant for dynamical electroweak symmetry breaking can have dif- 
ferent patterns of spontaneous breaking of the global symmetry than QCD. The 
resulting Goldstone Bosons, or pseudo-Goldstone bosons in the presence of mass 
terms, are thus in different manifolds and the low-energy EFT is also different. 

In this paper we only discuss the same cases as in [1] where the underlying strong 
interaction is vectorlike and all fermions have the same mass. Three main patterns 
of global symmetry show up. A thorough discussion tree level or lowest order (LO) 
is [9]. With n fermions^ in a complex representation the global symmetry group is 
SU{n)L X SU{n)ji and it is expected to be spontaneously broken to the diagonal 
subgroup SU{n)v- This is the direct extension of the QCD case. For n fermions 
in a real representation the global symmetry group is SU (2n) and it is expected to 
be spontaneously broken to SO{2n). In the case of two colours and n fermions in 
the fundamental (pseudo-real) representation the global symmetry group is again 
SU (2n) but here it is expected to be spontaneously broken to an Sp{2n) subgroup. 
Earlier references are [10, 11, 12]. Some earlier work for the complex case and the 
pseudo-real case at NLO can be found in [13, 14, 15]. 

In the remainder of this paper we refer to the complex representation case as 
complex or QCD, the real representation adjoint or real and the pseudo-real 

representation case as two-colour or pseudo-real. In [1] we extended the construction 
of the general Lagrangian to NLO^ including the divergence structure. The NNLO 
for the QCD case is in [16] and the divergence structure in [17]. The Lagrangian 
constructed in [16] is with the changes discussed in [1] and in Sect. 2 also a complete 
Lagrangian for the other two cases but we have not shown it to be minimal nor 
calculated the divergence structure. 

We do not repeat the discussion of the three different cases at the underlying 
fermion (quark) level. This can be found in [9] and [1], Sect. 2. In Sect. 2 we quote 
the structure of the effective field theories for the three cases but we again refer to [1] 

^We use n rather than Np for the number of flavours since it makes the formulas shorter. 
^References to some related work can be found in [5]. 



for more details. Sect. 3 discusses in detail the general structure of the amplitude. 
The amplitude can be expressed in terms of two functions B{s,t,u) and C{s,t,u) 
which are generahzations of the amplitude A{s,t,u) in TTTr-scattering [18]. We work 
out the possible intermediate states using the relevant group theory and using a 
projection operator formalism obtain the amplitudes in the different channels. The 
results for the amplitude arc discussed in Sect. 4 and for the scattering lengths in 
Sect. 5. Here we present some representative numerical results for the scattering 
lengths as well as some large n relations between the different cases. The lengthier 
formulas at two-loop order are given in an appendix. This work needed a few more 
integrals at intermediate stages than [19, 20], these are given in App. D. In Sect. 6 
we summarize our results. 



2 Effective Field Theory 

2.1 Generators 

The notation for the three cases can be brought in a very similar form. More details 
can be found in [1]. The Goldstone boson live on a manifold G/H where G is the full 
global symmetry group and H is the part that remains unbroken after spontaneous 
symmetry-breaking. We label the unbroken generators as and the broken ones as 

The space SU{n) x SU{n) / SU (n) is isomorphic to SU{n) so we use the X" as the 
generators of SU (n) for the QCD case. They are traceless, hermitian nxn matrices. 

The adjoint or real case has the generators in SU {2n) / S0{2n) where the broken 
generators satisfy 

JsX- = {X-f Js , with ^5 = ( 5 Q ) . (1) 

/ is the nxn unit matrix and the superscript T indicates the transpose. The X" 
are traceless, hermitian 2n x 2n matrices in this case. Multiplying (1) with Js from 
left and right leads immediately to 

X^Js = Js{X'^f. (2) 

The two-colour or pseudo-real case has the generators in SU{2n)/Sp{2n) where 
the broken generators satisfy 

JaX'' = (X^'f J A , with Ja=^] ) . (3) 

The X" are traceless, hermitian 2n x 2n matrices also in this case. Multiplying (3) 
with J A similar to above gives 

X'^Ja = J A {X'^Y ■ (4) 
The unbroken generators satisfy 

S0{2n): T''Js + JsT""^ = 0, 

Sp{2n) : T''Js + JsT"^ = . (5) 



This allows in both cases to derive using J = Js or J = Ja respectively: 

h)j^Jh^ with h^expik'^T^ . (6) 
We always use generators normalized to one: 

(T"T^) = = (7) 

(A) = tr p{A), is the trace over the flavour indices. This is over n for the QCD case 
and 2n for the real and pseudo-real case. 

During the course of the calculation, we often have to sum over the Goldstone 
Bosons. These sums can be easily performed using 

complex : 

= {A){B) --(AB) , 
n 

= (AB) --{A){B) . 

Th 

Real : 

{X'^AX'^B) = ^-{A){B) + ^{AJsB^Js) - ^{AB) , 

{X-A){X-B) = l{AB) + \{AJsB^Js)-^{A){B). 
Pseudoreal : 

{X'^AX-B) = \{A){B) + \{AJaB^Ja)-^{AB), 

= \{AB)-\{AJaB^Ja)-^{A){B). (8) 

There is a relation that the broken generators satisfy for the real and pseudo-real 
case. 

(x"x^ . . x^x') = . . x^'X") . (9) 

The proof for the real case is 

{x^xK-.x^x^) = (x"x^..x^x'J|) 

= {X^X^ ...X^JsX^'^Js) 

= (x"x^ . . j^x^^^x'^ J5) 

= {JsX'^^X'^ ...X'^^X'^'Js) 
= (X°^X''^...X^^X'^) 

= ((X'X^..X''X")^) 

= (X'X'^ . . . X^X") (10) 

The pseudo-real case is proven by replacing J| by —J\ and following the same steps. 
(9) is also the reason why the Lagrangian in [16] is not minimal for the real and 
pseudo-real case. 

In the group theory references there is a conjecture mentioned that to get from 
S0{2n) to Sp{2n) it is sufficient to take n —n. This feature is indeed visible in 
most of our formulas. 



2.2 Lagrangians 

As described in more detail in [1] we can write the Lagrangians in the three cases in 
a very similar way. The Goldstone Boson manifold G/H is parametrized by 

u^exp(^-^<j)Y (l> = rX\ (11) 

These transform under the symmetry transformation in the QCD case for Ql x 
qr e SU{n)L X SU{n)R as 

u gRuh{gL, Qr, 0)^ = h{gL, Qr, (t>)ugl . (12) 

h is the socalled compensator field and is defined by (12) and is also an SU (n) matrix. 
This can be derived from the standard general formulation [21] as done in [1]. For a 
transformation in the conserved part of the group we have that Ql = Qr = Qv and 
h = gv- 

The notation for the other two cases is directly that of [21]. A symmetry trans- 
formation g & G — SU (2n) transforms u as 

guh{g,(j)y , with = exp (i/i'^T") . (13) 

I.e. /i is in the unbroken part H of the group. In case the transformation g is in the 
conserved part of the group, g & H, we have that h — g. 
We can now define the quantities 

Under the group transformation in all cases we have Uf^ — > hu^h^ and can be used 
to define a covariant derivative. 

V^Wj. = dij,Ui, + F^Wj, - li^F^ hV^u„h) . (15) 

In [1] we also showed how the external fields can be included in a similar way as for 
the QCD case in [22, 13]. In particular the quark masses can be put in a quantity 

X± that transforms as x± ~^ hx±h^ ■ 

The lowest order Lagrangian takes on the standard form 

>CLo = ^(iiX + X+) (16) 
for all three cases and the same is true for the NLO Lagrangian. 

+L3{u^'uXuu) + L4(m%)(x+) + L5(m%x+) + Le{x+)^ 
1 

2' 

We have kept only the terms contributing to meson- meson scattering in (17). 

The NNLO Lagrangian is known for the complex or QCD case [16] as well as its 

divergence structure [17]. The same Lagrangian with the changes mentioned above 
is complete for the other two cases but probably not minimal. We have nonetheless 
chosen to leave the contributions from those terms in the results quoted here. 



+L,{x-r + -Ls{xl + x'-). (17) 



2.3 Renormalization 



We use the standard renormalization procedure in ChPT [22, 13] with the extension 
to NNLO described in great detail in [20, 17]. The divergences at NLO are canceled 
by the subtractions as calculated in [1]. At NNLO the divergences for the QCD case 
are canceled by the subtractions calculated in [17]. The other two cases satisfy all 
the expected constraints. Nonlocal divergences fully cancel, the e parts of the loop 
integrals as defined in App. D always cancel and the double divergences satisfy the 
Weinberg relations [17]. 

As usual in ChPT we apply the M^* scheme of dimensional regularization, in 
which the bare LECs Lj are defined as 

L, = M'^-^[r,A + LI(/.)] (18) 

Where the dimension d — A — 2e, and 

Inc = -i[ln47r + r'(l) + 1] . (20) 
2 

The coefficients Fj can be found in [22, 1] for the complex and in [1] for the real and 
pseudo-real5'p(2n) case. 

The NNLO terms can be made finite with the subtractions 



Ki = {cjj) 



2{d-4) 



(21) 



The coefficients rf \ F-^-* and F^^^^ for the complex case have been derived in [17]. 
For the real and pseudo-real case, the results do not exist. We have checked that all 
remaining divergences are local and can thus be subtracted. 



3 General results for the amplitudes 

3.1 TTTT case 

The TTTT scattering amplitude, which correspond to the QCD case with n = 2 is well 
known. Due to crossing and the possible SU{2) (isospin) invariants the amplitude 
can be written as [18, 23] 

M^^(s, t, u) = 5»^5"'^A(s, t, u) + S'^^S^Ait, u, s) + rt^"A(ii, s, t) . (22) 

The function A{s, t, u) is symmetric under the interchange of t and u. 

The possible states of two pions are isospin 0, 1 or 2. The amphtude for the three 
channels are given by [23] 



T\s,t,u) 
T\s,t,u) 



= 3A{s,t,u) + A{t,s,u) + A{u,t,s), 
= A{t, s,u) — A{u, s,t), 

= A{t,s,u) + A{u,s,t). (23) 



Where / is isospin, and Pi is the projection operator on isospin / They satisfy the 
relation 

M^,{s,t,u) = J2 T'{s,t,u)Pj. (24) 

7=0,2 

and 

T\s, t, u)Pi (no sum) = PiM^^{s, t, u) . (25) 

In the remainder of this section we will generalize these results. (22) is generalized 
in terms of two functions in Sect. 3.2. The possible intermediate states and the 
corresponding amplitudes are derived for the three cases separately in the last three 
subsections of this section. 



3.2 General amplitude 

The amplitude for meson-meson scattering is given by 

(0^(Pc)0'(Pd)|0"(Pa)0'(Pd)) = M(s,t,u) . 

The Mandelstam variables s, t, u are defined by 
These satisfy 



u 



(26) 



ipa-PdflMl^,,. (27) 



(28) 



We have chosen here to use the dimensionless versions in order to simplify later 
formulas. 

The flavour structure of the amplitude for meson-meson scattering can be de- 
scribed by constructing all possible invariants from the four corresponding generator 
matrices X^, e = a,b, c, d. Taking into account that (X^) = there are 9 invariants 
possible 



{X^X''X^X'^), 
(X'^X'^X^X^), 
{X''X^){X''X'^), 



{X'^X'^X^X''), {X^X'^X^X^), 
(X'^X^X'^X^), (X^X^X^X''), 
(X"X'=)(X''X'^), (X"X<^)(X''X^ 



(29) 



Under charge conjugation X" — > X"-^. This means that the amplitudes multiplying 
the first row in (29) must be the same as those multiplying the second row.^ As a 
result the full amplitude can be written in terms of two invariant amplitudes B{s, t, u) 
and C{s, t, u). 



M{s,t,u) 



Bis, t, u) 



+ 
+ 



(X^X^'X^X'^) + (X^X'^X'^X'') 

(X^X^X'^X^) + (X^X^X'^X^)] B{t,u,s) 
(X"X'^X^X^) + (X^X^X^X'^)] B{u,s,t) 



+S''''S'"^C{s, t, u) + S'^'S'^Cit, u, s) + S'^'S'^Ciu, s, t) 



s;acxbd 



adsbc/ 



(30) 



^Alternatively use (9) for the real and pseudo-real case. 



The flavour structure also implies that 



B{s, t, u) = B{u, t, s) C{s, t, u) = C{s, u, t) . 
For = 3 there is the Cayley-Hamilton relation 



{x^x^x'x'^) = {x''x^){x''x" 

6 permutations 3 permutations 



(31) 



(32) 



which allows for an ambiguity in the split of B and C . For n = 2 we can perform all 
the traces with four matrices X" in terms of Kronecker deltas. 

The relation between the general amplitudes and the vrvr scattering case (22) is 



A{s, t, u) = C{s, t, u) + B{s, t, u) + B{t, u, s) - B{u, s, t) 



(33) 



Note that the property (31) insures that A{s, t, u) is symmetric under the interchange 
of t and u as it should be. The form (22) holds also for any set of pions. I.e., taking 
any SU (2) subgroup of the unbroken group and any three of the pseudo Goldstone 
bosons that form a triplet under such a group, one can rewrite the general amplitude 
(30) into (22) using (33). 

3.3 QCD case: channels and amplitudes 

The Goldstone boson transform under the conserved part of the group, SU{n), as 

^ h(f)h^ . (34) 

This means that they are in the adjoint representation of SU {n). For the n = 2, 3 case 
we have an isospin triplet under SU{2) and an octet under SU{3). The intermediate 
states for these are well known: 



SU{2) : 3 ® 3 
SU{3): 8® 8 



1©3©5 (or / = 0,1,2), 
1 © 85 © 8a © 10 © To © 27 . 



(35) 



The group theory for SU {n) can be done in many ways. One is via Young diagrams 
and the second using tensor methods. We will do both. The SU (n) case derived here 
was in fact known [24] and our results are in agreement with his. Young diagrams 
for SU{n) are explained in [25] page 370. The Young diagrams for SU{n) give 

~n © rn = ■ © rn © rn © n~n © rrn © rn © i i i i i (36) 



□ 



□ 



Note that the : stand for n — 5 boxes. In terms of free indices the right hand side 
of (36) is no indices (singlet), twice one upper and one lower index (adjoint). The 
remaining four have all two lower and two upper indices, where the upper indices are 
produced from the columns with length n — 1 and n — 2 boxes using the Levi-Civita 
tensor e*i---*". For these we use the notation R\ where X = S,A indicate whether 



the lower indices are symmetric or antisymmetric and Y = S,A the same for the 
upper indices. The decomposition (36) can thus be written as 



Adj. (g) Adj. ^ Ri®Rs®Ra® +Rs^ ® Ra ® ® ^s^ (37) 

The order of the irreducible representation on the right-hand side is the same in (36) 
and (37). 

If we have a particular two-meson state 0"(pi)0''(p2) wc have to write it in terms 
of states that belong to the irreducible multiplets to obtain the amplitudes for the 
different channels. We use a simplified notation below with A = X", B = X^, C = X'^ 
and D — X'^ for many of the terms. All traces connecting a lower with an upper 
index must vanish. 

• Rj: singlet representation. All indices should be contracted, so this must be 
proportional to AjBj = (AB). Summing over the n^ — l states that are present 
tells us that the correct normalized state is 

Rj = -j^= J2{X''X'')r<p' . (38) 
Vri' - 1 a,b 

A projection operator p^'^'^'^ that projects on the singlet component is 

pabcd ^ _J_ ^^^^ ^^j^^ _ (39) 
Th -L 

It can be checked that this is a projection operator 

pabcdpcdef _ pabef {'iO) 

using (8) and Ea,6^a6 = - 1. 

• Rs'- adjoint symmetric representation. This needs in the end an upper and a 
lower index and must be traceless. We choose here to split up the two possible 
contractions of the indices of A and B in way that is symmetric under the 
interchange of A and B. 



n 



■s, 



^ Y 2(^2 - 4) 



AfBl, + B]^Al,--6]{AB) 
n 



(41) 



The last term is needed to make Rg traceless. The normalization can be worked 
out by checking the normalization of a particular state or by checking that the 
projection operator^ 

Ps^Rs{A,B)iRs{C,D)i (42) 

has the correct normalization, its square is equal to itself. This leads finally to 
the projection operator 



Ps 



{{AB + CD){CD + DC)) - -{AB){CD) 

Th 



(43) 



2(n2 - 4) 

Suppressing the indices we get similar to (40) P| = Ps and PsPi — PiPs — 0. 



*We suppress the superscript abed from now on. 



Ra'- adjoint anti-symmetric representation 

(RaTj = ^ {AfB\ - Bf^^) . (44) 

Ra{A, B) is traceless and antisymmetric in A and B. The projection operator 
corresponding to this is 

Pa = ^{{AB - BA){CD - DC)) . (45) 

Kg-, symmetric for lower indices and antisymmetric for upper. A state here 
corresponds to 



1 



n 



-61 {AlBl^-ATBl 



I m 



(46) 



-(i -H- j) {k ^ I) - {i ^ j,k ^ I) . 
The projection operator on this type of states is 

Ps. = (ii^(^,B))«(«|(C,i5));' 

= ^[{{AB-BA){CD-DC)) + n({ACBD) 

-{ADBC)) +n[{AC){BD) - {AD){CD))] . (47) 

R^: antisymmetric for lower indices and symmetric for upper. A state here 
corresponds to 



kl 



9 ^IBi - -SI {AlBr - A^bQ 

+ (^ j) - {k -fr^ I) - {i -fr^ j, k -fr^ I) . 

The projection operator on this type of states is 

Pas = —\((AB-BA)(CD-DC))-n({ACBD) 
An L ^ 

-{ADBC))+n({AC){BD) - {AD){CD))\ . 



(48) 



(49) 



Rg\ symmetric for both upper index and lower index. The states are 



{Ps)ki — 2 



:SlSi{AB) 



(n + l)(n + 2) 
+{i -H- j) + (/c -H- Z) + (i -H- j, /c -H- 1) 

and the projection operator is 

Pss = ^^^^{{AB + BA){CD + DC)) + \{{ACBD) 

+{ADBC)) + 1{{AC){BD) + {AD) {CD}] 
1 



(50) 



2(n + l)(n + 2) 



{AB){CD) 



(51) 



R^: antisymmetric for both upper index and lower index. The states are 

1 

n-2 

1 



(n- l)(n-2) 
-{i ^ j) - {k ^ I) + {i ^ j,k ^ I) (52) 
and the projection operator is 

-1 



Paa — 



-{{AB + BA){CD + DC)) - U{ACBD) 



4(n-2)'' " 4 

+{ADBC)) + U{AC){BD) + {AD) {CD)) 



4 

^ {AB){CD). (53) 



2(n- l)(n-2) 

The projection operators agree with those of [24] and using (9) can be shown to 
satisfy Pj-P^i — P^drr' for r the various representations. One last check is that 

Y,Pr^{AC){BD). (54) 

r 

The right-hand side is the unit operator when acting on the product of two states in 
the adjoint representation. It can also be seen that i?/, Rg^ Ra and i?| are symmetric 
under interchanging A and B while Ra-, Ra and Rg are antisymmetric. 

The amplitude in the different intermediate states can now be extracted from the 
general amplitude in two equivalent ways. We can pick a state Rr in a representation 
r and get it via 

Tr^{Rr\M{s,t,u)\Rr), (55) 

or apply the projection operators on the full amplitude with 

PrTr = PrM{s, t, u) . (56) 

Both methods give as expected the same result but the second one is much easier 
to apply. For the first method it is best to choose a state where the terms with 5 
functions are not present. E.g. for the four index representations take a state Rr 
with i — l,j — 2,k — 3,1 — A. For evaluating (56) one can use (8). 

Ti = 2(n--][B(s,t,u) + B(t,u,s)]--B(u,s,t) 
\ nj n 

+(n^ - l)C(s, t, u) + C{t, u, s) + C(u, s, t) , 

/ 4\ 4 
Ts = (n- - )[B{s,t,u) + B{t,u,s)]- -B{u,s,t) 

+C{t,u,s) + C{u,s,t) , 
Ta = n[—B{s,t,u) + B(t,u,s)] + C{t,u,s) — C{u,s,t), 
TsA = C(t,u,s) - C{u,s,t) , 
Tas = C(t,u,s) - C{u,s,t) , 
Tss = 2B{u,s,t) + Cit,u,s) + C{u,s,t) , 

Taa = -2B{u,s,t) + C{t,u,s) + C{u,s,t) . (57) 



They satisfy the relation similar to (24), 



M{s,t,u) = Y,Tr{s,t,u)Pr . (58) 

r 

One also notices that Tsa = Tas in general from (57). 



3.4 Real case: channels and amplitudes 

In this subsection we work out the possible two meson intermediate states for the 
case of SU{2n) / S0{2n) . One problem is that the mesons transform under S0{2n) 
as — )■ hcl)^^ . The matrices h G S0{2n) in the embedding introduced here do 
not simply satisfy hh^ = 1 either. In our case the S0{2n) is instead defined as 
hJsh^ = Js- The easiest way to obtain objects that appear in the usual way is to 
note that using (6) 

(j)Js -> h(j)hKjs = HJsh^ (59) 
and that an invariant trace on these object needs an extra factor of Js- E.g. 



(0V5) J5(0V5) ^ KrJsW JsK<p'JsW = K<p''Js)Js{<P'JsW ■ (60) 

Keeping that in mind we can use the standard way of dealing with S0{2n). Note that 
{(f)Js)^ = Js4>'^ = 4>Js so the Goldstone bosons live in the symmetric representation 
of S0{2n). 

The method of Young tableaux has been generalized to S0{2n) [26]. Putting 
together two symmetric representations gives 



(61) 



We can write this in the form 



Sym. ®Sym. = Ri® Ra® Rs® Rfs © Raia © Rms ■ (62) 

The states are given in Tab. 1. They are made traceless but remember that indices 
of (f)°'Js and (p^Js are always contracted with J5. 

As an example the singlet representation Rj is proportional to 

{rJs)^JJs^kJsA(p'Js)kl = (^0') • 

The two-index antisymmetric representation Ra is, now using A, B for 0", (^*, 

{AJsU{BJs),iJski - {AJs)jk{BJs)uJski = {ABJ - BAJ),j (63) 

where we heavily used J| = 1 and the fact that AJs and BJs are symmetric matrices. 
One can also easily check that the trace {ABJs — BAJs)ijJsij vanishes. 

Then Ra and Rs are antisymmetric respectively symmetric both under i ^ j 
and A ^ B. The remaining ones are the four-index representations. Rps is fully 
symmetric in all indices and under A ^ B. The two remaining representations have 
a mixed symmetry in the indices but are antisymmetric respectively symmetric under 
A^ B. 



Ri 


y/ (2n-l)(n+l) 


Ra 


{ABJ - BAJ)ij 


Rs 


{ABJ + BAJ)ij - ^Jij{AB) 


Rps 


{AJ),j {BJ),i - -±-^[j^^{ABJ + BAJ),i] + ^^-^^J^^J,i{AB) 
+ {ijkl -H- ikjl) + {ijkl O iljk) + {ijkl -H- klij) 
+ {ijkl ^ jlik) + {ijkl ^ jkil) 


Rma 


{AJ)ij {BJ)ki - {AJ)ki {BJ),j -^[j,,{ABJ - BAJ)ji 
+Jjk{ABJ - BAJ)u + Ju{ABJ - BAJ)jk + Jji{ABJ - BAJ)ik\ 


Rms 


{AJ)ij {BJ)m + {AJ)ki {BJ)ij - {AJ)ik {BJ)ji - {AJ)ji {BJ)ik 
+2i^ [jiMBJ + BAJ)ki + Jki{ABJ + BAJ)ij 
-Jik(ABJ + BAJ)ji - Jji{ABJ + BAJ)ik] 

(n-l)(2n-l) {jtjJkl JikJjl) {^B) 



Table 1: The intermediate states for the real or adjoint case, SU{2n) / S0{2n). The 
notations A, B stands for 0" and (f)^ and J is Js everywhere. 



Pi 




Pa 


-2i^)UB-BA){CD-DC)) 


Ps 


2(n-i)(nH-2)(((-4^ + BA){CD + DC)) - ^{AB){CD)) 


Pes 


l[(^^+lkw){^B){CD) + {AC){BD) + {AD){BC) 
+2{ACBD + ADBC) - ^{{AB + BA){CD + DC))] 


Pma 


^^{{AB - BA){CD - DC)) + \{{AC){BD) - {AD){BC)) 


Pms 


l[in^mn-iMB){CD) + 2({AC){BD) + {AD){BC)) 
-2{ADBC + ACBD) - ^{{AB + BA){CD + DC))] 



Table 2: The projection operator for the different intermediate states for the real or 
adjoint case, SU{2n) / S0{2n). 



The corresponding projection operators are given in Tab. 2. These can be ob- 
tained by contracting the indices with J5 of the states once with A, B and once with 
C,D. 

We can now use Tabs. 1 and 2 to project out the amphtudes in the various 
channels using (55) or (56). The results are 

Ti = ^{2n-l){n + l)[B{s,t,u) + B{t,u,s)] + ^{n-l)B{u,s,t) 

+ {2n -l){n + l)C{s, t, u) + C{t, u, s) + C{u, s, t) , 
Ta = -il + n)[B{s,t,u)- B{t,u,s)] + C{t,u,s)-C{u,s,t), 

Ts = -{n-l){n + 2)\B{s,t,u) + B(t,u,s)] + -{n-2)B(u,s,t) 
n n 

+C{t,u,s) + C{u,s,t) , 
Tfs = 2B{u,s,t) + Cit,u,s) + C{u,s,t), 
Tma = C{t,u,s) - C{u,s,t) , 

Tms = -B{u,s,t) + C{t,u,s) + Ciu,s,t). (64) 



3.5 Pseudo-real case: channels and amplitudes 

In this subsection we work out the possible two meson intermediate states for the 
case of SU{2n)/ Sp{2n). One problem is that the mesons transform under Sp{2n) as 
in — 7- h(f)h'' . The easiest way to obtain objects that appear in the usual way is to 
note that using (6) 

4>Ja h(j)h^JA = h(j)JAh^ (65) 

and that an invariant trace on these object needs an extra factor of Ja- Keep- 
ing that in mind we can use the standard way of dealing with Sp{2n). Note that 
{(PJa)"^ = Ja4>^ = —Ja4>^ = <pJs so the Goldstone bosons live in the antisymmetric 
representation of Sp{2n). 

The method of Young tableaux also is developed for Sp{2n) [27] . Putting together 
two antisymmetric representations gives 



(66) 



We can write this in the form 



Asym. ®Asym. = Ri® Ra® Rs® Rfa ® Rma ® Rms ■ (67) 

The states are given in Tab. 1. They are made traceless but remember that indices 
of (P^'Ja and (P'^Ja are always contracted with Ja- 

The representations are the singlet representation, symmetric under A -H- 5, Ra 
which is antisymmetric under the interchange i ^ j but symmetric under A ^ B 
and Rs which is symmetric under the interchange i ^ j but antisymmetric under 
A ^ B Let us show the latter on Ra The two-index antisymmetric representation 



Ri 


^(2„+l)(„_l) \ / 


Ra 




Rs 


{ABJ - BAJ)ij 


J^FA 


{AJ)ij {BJ)ki + ^^JiA^BJ + BAJ\i - ^^-^J,^J,^{AB) 
-{ijkl <H- ikjl) + {ijkl -H- iljk) + {ijkl <H- klij) 
+{ijkl -H- jlik) + {ijkl -H- jkil) 


Rma 


{AJ)ij {BJ)okl - iAJ)ki {BJ)oij - 
+Jji{ABJ - BAJ)ik - MABJ - BAJ)jk - 


'jikiABJ - BAJ)ji 
- JjkiABJ - BAJ)a] 


Rms 


{{AJh {BJ)ui + (AJ),fc {BJ),i - 

-Jik{ABJ + BAJ)ji\ + {ij < 

' 2(n+l)(2n-l) ("^^i"^*^' '^^kJjl 


J^j{ABJ + BAJ)ki 

kl)} 
){AB) 



Table 3: The intermediate states for the pseudo-real or two-colour case, 
SU{2n)/ Sp{2n). J means J^. 



Ra is, now using A,B for 0", <^'', 

{AJA)ik{BJA)jlJAkl — {AJA)jk{BJA)ilJAkl 

— —{AJA)ikJAkl{BJA)ljJAkl — {B J A)ilJ Alk{AJ A)ki 

= {ABJ + BAJ)ij, (68) 

where we used Ja — —1 and the fact that A J a, BJa and J a are antisymmetric 
matrices. 

The remaining ones are the four-index representations. Rfa is fully antisymmet- 
ric in all indices and symmetric under A ^ B. The two remaining representations 
have a mixed symmetry in the indices but are antisymmetric respectively symmetric 
under A-<r^ B. The states are give in Tab. 3. 

The projection operators can be constructed by contracting all indices with Ja 
of the states once with A, B and once with C, D. The results are given in Tab. 4. 

We can now use again (55) or (56) to obtain the amplitudes in the different 
channels. The results are very similar to the real case and read 

Tf = -{2n+l){n-l)\B{s,t,u) + B{t,u,s)\--{n + l)B(u,s,t) 
n n 

+(2n + l)(n - l)C(s, t, u) + C{t, u, s) + C{u, s, t) , 

Ta = -{n + l){n-2)[B{s,t,u) + B(t,u,s)]--(n + 2)B(u,s,t) 
n n 

+C{t,u,s) + C{u,s,t) , 
Ts = {l-n)[B{s,t,u)-B{t,u,s)] + C{t,u,s)-C{u,s,t), 
TpA = -2B{u,s,t) + C{t,u,s) + C{u,s,t), 
Tma = C{t,u,s) - C{u,s,t) . 
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Pa 


2(n+i)(n^2)(((^^ + BA){CD + DC)) - '-{AB){CD)) 


Ps 


-^{{AB-BA){CD-DC)) 


PpA 


1 

6 

-2( 


'^-^{AB){CD) + {AC){BD) + {AD){BC) 
ACBD + ADBC) - ^{{AB + BA){CD + DC))] 


Pma 


^{{AB - BA){CD - DC)) + \{{AC){BD) - {AD){BC)) 


Pais 


1 

3 

+{A 


^—^{AB){CD) + {AC){BD) + {AD){BC) 
CBD + ADBC) - ^^^^{iAB + BA){CD + DC))] 



Table 4: The projection operator for the different channels for the pseudo-real or 
two-colour case, SU{2n)/ Sp{2n). 




(1) (2) (3) (4) 

Figure 1: The leading order and next-to leading order for meson- meson scattering 
00 — )■ 00. The filled circle is a vertex from C2, and the filled square is a vertex from 
£4. 

Tms = B{u,s,t) + C{t,u,s) + C{u,s,t) , (69) 



4 Results for the amplitude M(5, u) 

We have rewritten the amplitudes here in terms of the physical decay constant -Fphys 
and mass M^^iys- The relation of these to the lowest order results can be found in [1]. 
We also use the abbreviations 

X2 = — , L = -m — , TTifi = ^7. 70 

Fp\y, 167r2 /i2 I67r2 ^ ' 

In addition we have often used (28) to simplify the expressions. 



4.1 Lowest order 

The lowest order result comes from the simple tree- level diagram (1) of Fig. 1. 



BLois,t,u) = X2[--t+lj , CLo{s,t,u)^0, (71) 

for all cases. This reproduces using the relation (33) Weinberg's result [18] for tttt 
scattering 

Alo = X2 (s - 1) . (72) 

4.2 Next-to- leading order 

The next-to-leading order contains the 3 diagrams (2-4) in Fig. 1 in addition to 
wave-function-renormalization. The divergences from loop diagram (3) and (4) can 
be canceled by the bare low energy constants (LECs) of £4 in the diagram (2). 

The functions B{s, t, u) and C(s, t, u) can be calculated from the one-loop graphs 
shown in Fig. 1(2-4) and wave-function renormalization. It also contains terms from 
rewriting the lowest-order result in the physical mass and decay constant. 

The functions B{s, t, u) and C(s, t, u) can be rewritten in the form 

B{s, t, u) = xl [Bp{s, t, u) + Bs{s, t-u) + Bs{u, t - s) + Brit)] , 
C{s,t,u) = xl[Cp{s,t,u) + Csis) + CTit)+CTiu)] . (73) 

Bp{s,t,u) and Cp{s,t,u) are the polynomial part, the remaining pieces are often 
called the unitarity correction. This can be proven using an extension of the methods 
of [28]. 

Using (28) wc rewrite the polynomial part in its simplest form satisfying the 
symmetry constraints (31): 

Bp{s,t,u) — ai + + a^t"^ + a4{s — , 

Cp(s,t,u) = iSi + P^s + Pss^ + P^it-u)^ (74) 

The polynomial part for the three cases is give in Tab. 5. 

The unitarity correction is given in Tab. 6. We noticed that the C functions for 
the S0{2n) and Sp{2n) case are the same. 

4.3 Next-to-next-to-leading order 

There are 13 diagrams at next-to- next-to-leading order shown in Fig. 2. We have 
checked that the nonlocal divergence cancels for all three cases and that for the 
complex or QCD case the result is fully finite with the subtractions calculated in 
[17]. 

The diagrams (15) and (16) are often called the sunset and vertex or fish diagram 
respectively. These require the most difficult integrals. At intermediate stages we 
needed more integrals than those calculated for [19, 20]. They were calculated with 
the methods of [29] and are given in App. D. 



B{s,t,u) = xl[Bp{s,t,u) + Bs{s,t-u) + Bs{u,t-s)+BT{t)] , 
C{s,t,u) = xl[Cp{s,t,u) + Cs{s) + CT{t)+CT{u)] . (75) 



QCD: SUin) x SU{n)/SU{n) 


Bp{s,t,u) 


ai 


f L + f TTie + 16L^ + - \nL - fnyrie 




Ol2 






«3 


^3 + 4i^5 - ^nL - ^riTTie 






^3 - - ak^^ie 


Cp{s,t,u) 




32(^-L^ + L^)-^(L + 7ri6) 
16L^ - 32 L\ 






-\L + 2LI + SL^ - fTTie 




Pa 


2^2 - 1^ - I^Tie 


Adjoint: SU{2n) / S0{2n) 


Bp{s,t,u) 


ai 


16(i^0 + ^8)-(i + f + f -^)-16 






-16L?S - 4L'5 + (ft + 1) + (1 + 1^) vrie 






4i^o + i^S-(| + f|)^-(i + tl) ^16 




0:4 


^3 + (2! ~ H) ^ + (lli ~ §) 


Cp{s,t,u) 




32(Li-L^ + L[;)- 2^(L + 7ri6) 
16 {LI - 2L\) 




As 


8L^ + 2L2 — ^TTie — ^1/ 


Two-colour: SU{2N)/ Sp{2N) 


Bp{s,t,u) 




16(LS + LD + (i + i-f)L+(i|-f + i)7r,e 






-I6L5 - 4LE + (ft - 1) ^ + (1 - 1^) ^16 




«3 


4^S + i^3+(|-T|)^+(4-ti)^i6 




^4 


-^3 - (^ + 1) - (lii + ^16 


Cp{s,t,u) 


/34 


32(L'[-L^ + L^6)-2i.(L + 7ri6) 
16(L^ - 2U^ 
8LI + 2L2 — ^TTie — 

2-^2 - 16-^ ~ 16^16 



Table 5: The next-to-leading results for all three cases for the polynomial part. The 
coefficients are defined in (74). 




Figure 2: The next-to-next-to leading order diagrams for tttt — ?► tttt. The filled circle 
a vertex from £2, The filled square is a vertex from £4, and the open square is a 
vertex from C^. 



QCD: SU{n) x SU{n)/SU{n) 


Bsis,t-u) 
Brit) 


ns) 




f) {t - u)] 


Cs{s) 
Crit) 


Jis) + Is^) 
\J{t) {t-2f 




Adjoint: SU{2n)/SO{2n) 


Bsis,t- u) 
Brit) 


J{s) 


-2^ + 1 + 1^(^-1)^'+ (l-l) (t-u) 

\jm-2)^ 




Cs{s) 
Crit) 


As) + Is^] 
iJ(t)(t-2)2 




Two-colour: SU{2N)/ Sp{2N) 


Bsis,t-u) 
Brit) 


J{s) 


'-h-\' + U^ + ^)'' + U^ 

-|J(t)(t-2)2 


-1) (l-l) {t-u) 




Csis) 
Crit) 


As) + Is^] 
iJ(t)(t-2)2 





Table 6: The next-to-leading results for all three cases for the unitarity correction. 

The polynomial parts we rewrite using (28) in their simplest form satisfying the 
symmetry constraints (31): 

Bp{s,t,u) = 7i + 72^ + 73^^ + 74(s - -w)^ + 75^^ + 76^(s - ^i)^ 
Cp{s,t,u) = Si + S2S + Sss"^ + Siit - u)^ + 858^ + SQs{t - u)^ . (76) 

The coefficients in these polynomials as well as the functions in (75) are given in 
App. A. The FORM expressions can be downloaded from [30]. We stress once more 
that the result is fully analytical and expressed in terms of L and J. 

5 Scattering lengths 

The threshold parameters of general meson meson scattering are defined similar to 
those of nn scattering. First we calculate the amplitudes for the different channels 
using (57), (64) and (69) for each of the possible intermediate states of representation 
or channels /. 

The scattering amplitude for each channel / can be projected out using the partial 
wave expansion 

r/(s) = ^ j'^ d{cos9)Pe{cose)Ti{s, t, u) . (77) 

Near the threshold s = 4, we can expand the amplitude above the threshold using 
s — 4(1 -|- q^/M^) in the small three-momentum q. 

ReT^is) = q''[al + q%l + 0{q')], (78) 



where a| is the scattering length, and is the slope. 

In App. C, we give the expressions of the lowest partial wave scattering length for 
each channel in all three cases. As mentioned in Sect. 3, some channels are symmetric 
under A S, hence the lowest order partial wave is ^ = 0. The other channels are 
antisymmetric under A B, so that the lowest order partial wave is £ — 1. A and 
B are the incoming mesons here using the notation of Sect. 3. 

For the purpose of illustration, we plot the scattering length for the singlet and the 
fully-symmetric (fully-antisymmetric) channels as a function of the physical meson 
mass Mpj^yg. Since currently we do not have knowledge for the values of the low 
energy constants for these, we take the values of the of fit 10 of [31] for the 
complex or QCD case and half that for the other two as suggested by the large n 
relations discussed below. The values of the NNLO LECs we simply put to zero. 
We also choose the subtraction scale fx — 0.77 GeV and the physical decay constant 
Fphys = 0.0924 GeV. 

The singlet case for the complex case is shown in Fig. 3. We have divided the 
scattering length by n to make the lowest order similar for all cases. Plotted are 
n = 2, . . . , 5. One can see that for a given M^^^^ the convergence gets progressively 
worse for larger n. For n — 2 this corresponds to Uq. 

The fully symmetric case case for the complex case is shown in Fig. 4. Plotted are 
n — 2, ... ,5. One can see that for a given M^^iyg the convergence gets progressively 
worse for larger n. For n = 2 this corresponds to a^. The lowest order is independent 
of n. The NLO order is only mildly dependent on n while the NNLO part grows fast 
with n. 



5.1 Large n behaviour 

Looking at the lowest order expressions in App. C we notice immediately that the 
large n behaviour for the scattering lengths fall in three classes. 

The scattering length is of order n for the singlet and symmetric and asymmetric 
representation and what is more they are clearly related for the three cases: 

^0 1 complex 
"■0 [complex 
f^l [complex 

The symbol =lo means equality at lowest order. Do the relations (79-81) remain 
valid at higher orders? If we choose -Fp^ys oc n and Mpj^y^ independent of n, we find 
that it is indeed the case for (80,81). For (79) it is true, provided we set the NLO 
LECs LJ" of the real and pseudoreal case to half those of the complex case and the 
NNLO coefficients to 1/4 the complex case. The contributions are nonzero at the 
three orders for all of them. The subleading orders in 1/n are different. 

A second class is those which are order 1 in the coefficient of 2:2/77 at lowest order. 
For these we find 

_ ^FS\ _ 9„M5| _ ^ Zi 

w-Q [complex — "'0 |real — ■^"'0 Ipseudoreal — LO ^ ) v" / 



/, _ /, _ X2n 

C^olreal I pseudoreal LO „ ' \''^) 

71 O 

S\ -A I ^ /oA\ 

C^O Ireal — ^0 Ipseudoreal —LO T77 ; \°^) 

n io 

C^l Ireal (^l Ipseudoreal LO 777 ' \^^) 

TT 40 



n = 2 



n =3 




Figure 4: Scattering length Og"^ for the complex or QCD case, SU{n) x SU{n) / SU{n). 



n = 2 



n = 4 




Figure 7: Scattering length of Og/n for the pseudoreal or two-colour case, 
SU{2n)/Sp{2n). 



n = 2 



n = 4 




Figure 8: Scattering length 2aQ for the pseudoreal or two-colour case, 
SU{2n)/ Sp{2n). The factor of 2 included is because of the large n relation (82). 



"'Q Icomplex |real Ipseudoreal LO • v"'^/ 

TT iO 

We do indeed find that the relations are also satisfied at NLO and NNLO. In fact 
none of the scattering lengths in (82) has a leading n NLO correction to the lowest 
order result. This can be clearly seen in Fig. 4 where the NLO result is very similar 
in all plots. 

The third class is the amplitudes that vanish at lowest order 

„SA| _ AS\ _ n MA\ _ n„MA\ _ r, (oa\ 

^1 Icomplex ^1 Icomplex -^tl-^ |real ^Ll-^ Ipseudoreal LO " • v"^/ 

These are always suppressed by two powers of n compared to the first set of scattering 
lengths also at NLO and NNLO. The relations (84) are satisfied at NLO with same 
identifications of the LECs as above and almost at NNLO. The only terms that do 
not satisfy the relation are proportional to L^^L^. 

By comparing the plots shown one sees that the large n relations do predict the 
general behaviour but for n = 2 and n = 4 are not that accurate. 



6 Conclusions 

In this work we have presented the calculation of general meson-meson scattering for 
n flavours in a complex, real or pseudoreal representation of a strongly interaction 
gauge group. These are also referred to as QCD, Adjoint QCD and Two-color QCD 
and have as symmetry breaking patterns SU{n) x SU{n)/ SU{n), SU{2n) / S0{2n) 
and SU{2n)/Sp{2n) 

We first reviewed the effective field theories of these three different cases. Those 
theories can be written in a very similar form as discussed earlier [1]. We have 
extended the methods used for tttt scattering in ChPT [20] to all the present cases. 
At intermediate stages some more integrals showed up, we have calculated them and 
they are tabulated in an appendix. 

The amplitude can in general be written in terms of two invariant amplitudes 
which we called B[s,t,u) and C{s,t,u). These amplitudes can be written in terms 
of simpler functions and we have given their fully analytical expressions to NNLO. 



Since the long term use of our work is the study of scattering on the lattice for 
these alternative theories we have discussed the group theory involved and all the 
possible intermediate channels. We have derived the amplitudes in all these channel 
as a function of the invariant B and C functions. 

The expressions for the different channels we have not shown explicitly but we 
included expressions for scattering length of the lowest partial wave in all channels. 
We presented a few representative numerical results for the scattering lengths and 
discussed a series of relations between the different theories in the hmit of a large 
number of flavours n. 
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A Next-to-next-to leading order result 
A.l Complex or QCD 
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3 1 1 1 
~4^ ^ n2 ^ 24 I 



+m[ - j^^Lnh' - — + Ly + ^ + ^ + ^ + _w 

2rr 3 3 17 , 3 1 3 29tTi6S^ ^ r 2 2 

3 ° 3 ^ 576 288 192 144 
31Ls2 8L5s2 SL^s^ 28^5^ 13 

+^6 1 i ^ + 2Lls^-ALls^-Lls^-—Lns^ 

4 2 14 2 , rr- 2 35 2 2 , 19 3 557ri6S^ 
--LqUS - —L.^ns +L^ns - —n mes + ^^ttigS 

+ ^Ll{t - u)s^ + ^Ln(t - u)s^ - ^L^nit - u)s^ + '^Lln{t - u)s^ 

1 2 / X 2 47n7ri6(t-M)s2 437ri6(t - m^s^ 2QUr,s'^ 
n 7riQ[t — u)s 



432 ' 1728 1728 3n 

20L5s2 5 2 97Ls 4L5s l%L\s 32L^s 



— Ln's - + + + - 8L^s + 2L^s 



3n 36 72 

.7-r- 17Lns SLlns l&Llns ^ ^„ 
+16L^s + ALls - + + ^ ALlns + 8L^ns 

5 2 , riTiiQS ttiqS -kiqS ISvrieS 1^2/, n L{t - u)s 

+ 12" "^«^ + ^-^+2;^ + ^ + 48^" 

AUo{t-u)s 8L[{t-u)s Amt-u)s 2Ll{t - u)s ALl{t - u)s 

^ 3 3 ^ 3 3 3 

5(^^ u)s ^ _ ^-^g _|_ _]^r^^^ _ ^"^g _ -i,'^n(t — u)s — -Lln(t — u)s 

n . , , 253 , , IbbniQit - u)s Ls ULls 

3n 4n2 18 3 3 ^ ^ 18 6 

47ri6 , STTie , 27ri6 1 r 2/^ n L{t - u) l%Ll{t - u) 

-— + ^ + — - 12^" <*-")-— 6— + 3 

4La^ _ Ln(i - u) 4L;„(f - u) _ 11 _ 

3 12 3 48 ^ ^ 

59n7ri6(t ~ ''^) _|_ T^i^it — u) iriQit ~ u) ^ 37ri6(t — u) ^5L 80Lq 



144 8n 8n2 8 2n 3n 

80LS 16LS 48L'« 3L' 



J 5nt3 llt^ nt^ 3t^ nt t t 31t n 1 
^ ^^'[ 288" ~ 288 ^I6^T6~72^4^~4^~72'~12~2^ 

^2n2 ^ 12 I 



-, Ji9Lt3 sm^ i6m^ 2m^ 5 . 17 , 



96 18 

11, 2 1 2 SSTTiei^ 15m 64L^i 176L^i 
-36^"* + + _ + _ - 24L;t - ^ ^ 

- i6L;t - mt - 32m - t.Lit + + "^la^ + 



36 72 n n' 

27.,t ii + 4_68L^ 64L;^128L£^16I|^ 



2 



8 2n 2n2 18 ° 3 3 3 ^ 

7Ln 5n7ri6 27ri6 _|_ Svrie _|_ SSttis . L L 
n 'n? 18 n 



+16L^_l^_l^_^ + ^ + !^ + -_ii (90) 



, , , [ 3ns3 s^ 3s2 1 3 1 

r 13ns3 s3 23ns2 ^ ^ 1 ^ ^1 

+ '^^^\^88~ ^ 288 ~ ^iT ~72^T2"~2;^^4^2;?~;^^4j 

jLs^ 4LSs3 SL^s^ 5 . 3 8^, 3 85 3 

+ — |- + — |- - Y^^^s + 8L>s2 + -L^ns^ + ^^TTies' 

-l^-^-^ + 8Lis^ + ^-^ + 2L^,s^ + l^W 
288 72 3 ^3 3 ^72 

-32L>s^ - + l^Llns' - -nir^.s' + ^ _ _ 

3 12 16 n 3n 

-8L^s + 16L^s - + 32Lins + ^^^^^^ - 32L>s + 32L^ns 

18 3 

^35n7ri6S 77ri6S ^ SQttiss _^ Ls ^ 32L^s ^ I6L2S I6L4S 64LqS 
18 4n 72 n n 3n n 'i'n? 

''^^^ - + 32LI + ^ - 32L^ + 32L^ ■ '"^^ ''"^^ ■ 



3n2 3 ^ 3 ^ ^ n2 2n3 ^ 9 

L 2,2L\ 32L^ ^ 32L^ 2,2Uq 2L ^ SOU^ ^ 8OL5 16Lg 



n n 3n n n 3n2 3n'^ -n? 

i^Ll 7L 

H 

2^3 



48^8 7L 1 

+^ + 7r:^ (91) 



+J(t)l ^ + —Lnt^ nniat^ h UL^f 

^ ^{18 ° 3 144 288 288 18 ° 

H ^ + 2L^i2 ^ ^2 ^ ^ 24m ^ 

3 ^ 36 18 12 36 ° 3 

ATT, orr, , 31Lni 29n7ri6t , niat ISniat Lt 11 L , ^^^^ 
-4L5t-8V + ^^ + ^^ + -- + ^^---^ + 16Lo 

+l^ + 16L^-^-^-^-^ + -| (92) 
3 ^9 9 n6 nj 

A. 3 Pseudo-real or two-colour 

^-^^'^-^) = ^^^^H^^^^^"6 + 2;^+2;?r^-") 

^^^[ 48 18 48 ^ ' 36 24 24 

{t — u)n ^ s s{t — u) 7{t — u) 1 s (t — u) 



72 9 48 18 48 72 36 12n 12n 

1 s ^{t-u) 2 \ 



2n 12n^ 12n^ 3n^ j 

, , , (n^s^ ns^ 3s^ 3ns^ 1 9/ . ■? 1 / . ■? 

+ -64 + 64 - ^ - 576" <* " + 288"" " "^^ 

(t — u)s^ 9s^ 1 n, , n(t — u)s (t — u)s 3s n^it — u) 



576 32 72 ' ' 36 72 16 36 

nit-u) it-u) 1 1 1 



18 36 2n 2n2j 

, , , fn^s^ ns^ \ls^ VJrt^s^ llns^ 1 9/ x 9 

+A;i(s)<^ ^ \ \ nHt - u)s^ 

^ ^ \ 576 96 576 288 288 576 ^ ' 

1 , , , 5(t — ■u)s^ 3s^ n^s ns 1 , 

H nit - u)s^ ^ \ nHt - u)s 

144 ^ ' 576 32 12 144 96 ^ ' 

n(t — u)s (t — u)s s s 31s n n^it — u) 



32 48 8n 8^2 144 6 48 
n{t -u) 5(t-u) 3 1 1 1 
^ 24 ^ 48 ^ 4^ ^ ^ ^ 24 J 

J 5.3 19Ls3 3 ALls^ SU^s^ Lls^ 1 ^ 3 

+J(s)< Ln^s^ L^s^ H ^ \ Lns^ 

^ '[ 144 144 ° 3 3 3 96 

2rr- 3 3 17 2 3 1 3 297ri6S^ 7^22 



3 " 3 576 288 192 144 

+ ^ _ ^ _ ^ _ 2Lrs^ - AU^s^ + Us^ + i^Ln.^ 



4,r 2 14 2 , 2 35 2 2 19 2 557ri65^ 
--L^ns - —L^ns + L^ns - —n ttigs - ^nTn^s 

+^L(t - u)s' + hl,(t - u)s' + hl{t - u)s' - h^,(t - u)s' - hl(t - u)s' 



^ 1728 1728 3ri 3ri 36 ^ ~ 72 

4LSs IQLls 32^8 17Lns 8L[;ns 
^ + + + 8Lls - 2Lls + 16Lls - ALls + + ^— 

- '^^^^ + + 12^ ^^^^ - ^ + ^+ 2;^ + ^ 

1 , , L(t-u)s ALUt-u)s 8LUt-u)s AUJi - u)s 
+—Ln^(t - u)s ^ — ^ ^ i^- ^ H ^— 

48 ^ ^ 24 3 3 3 

2Lat-u)s 4:Ll{t-u)s LUt-u)s 1 , , 4 , , 
+^^^ ^ + 3 ^ + - u)s + -LSn(i - u)s 

2 1 17 253 

— Llnit — m)s Llnit — u)s H v?-Kie,{t — u)s nnic,(t — u)s 

33V 7351 ^ 216 ^ ' 864 ' 

155nie{t -u)s Ls QAL^s UL^s Ls 19L 16L[; 32L5 
^ 864 4^ ^ 3n ^ ~3n 4^ ^ 3 3~ ^ ^ 

5Ln SnTTie , 47ri6 , Sttig , 2ni6 1 r 2/. n L{t - u) 
-16i^8 + ^--^ + ^ + ^ + -9--^i^- (^--)-— ^ 

16L^(i-'u) 4L^(i-'u) Lnit-u) AUnit-u) U o , 

H ^ H ^ H ^ n^TTieit - u) 

3 3 12 3 48 ^ ^ 

59n7iiQ{t — u) 7iiQ{t — u) 7ri6(t — m) SttiqII — u) 5L 80Lq 8OL3 

^ 144 8n Sr^ ^ 8 2^ ~ 3^ 

^16L| _ 48L| _ 3L1 3 

n n ] 

B (t) - k (t)i ^^'^ t t t n 1 1 5 1 

T{t) - 3(t)| ^^ + ^ + y + ^ + ^ + Y8~T8~3n~3n2~18j 

, , J 3t3 9t 3 1 

+k2(t)l + + -} 

^ '[ 32 16 8 4J 

f 5nt3 11^3 3t2 nt t t 31t n 

^ '^[ 288 ~ 288 'l6'^T6'^72~4^~4^~72^12 

1 1 5 ) 

^2^^ 2^^ 12] 

J i9Lt3 8L':t=^ 2m^ 5,0 17 



3 



96 18 "3 

11, 2 1 2 557ri6t2 15m 64L':t 176L^t 16L^i 

H Lnt^ nvTigt^ + — + + 24L^t ^ ^ + ^- 

36 18 36 36 ° 3 3 3 

-16L^t + ALlt - 32Llt + 8L^t - ^ - _ ^^^^ _ Z^^at _ 27pt 



2n 2n2 18 " 3 3 3 ^ ^ 9 

BnTTiQ ^ 2-KiQ 2-KiQ 557ri6 L L 
n 'n? 18 n 



5ns s 7s 1 2 1 



5ns s s 1 111 

ZyoS"^ 5 O 8 Q 85 Q 

3 18 ^ 3 ^ 288 

^-8L^.^-^-^ + 2L^.^ + |w 

^ifir^' 2 7 2 SttisS^ 8L^52 8L-2s2 
f 16L,ns - -nvr^a. - ^ 



18 3 18 

Tttigs _ 897ri6S ^ Ls ^ 32L{s ^ 16L^s _ 16L^s _ 64LSs _ 64X^5 ^ L _ ^^^r 

4n 72 n n 3n n 3n^ 3n^ 3 ^ 



2n^ 9 n n 3n 
32K 32L^ 2L 80L:; 80L^ 16L^ 48L^ 7L 1 
+ - + — + ^ + ^ ^ + ^ + ^ 95 



^ , , , , Jnt'^ nt t t n 1 1 1 

^-W = '='W{^-S-ir-fe + 9+ 18 + 3^-18} 
, , J 3t^ 3t 1 1 

f 5nt=^ nt^ t2 ^ lit n 1 1 1 

^ ^^^'\ ~ 288 ~288^T6^T6~72^4^~72"~T2~2^^12J 

11 . 1 2 115*2 37£; 

— Lnt H nTrifii 

36 18 12 36 



IQLlt SlLnt 29nniQt tiiqI I'iniQt Lt 
— 41vc;t — oLiot H h r — — z — 

3 ^ ^ 36 72 n 72 2n 



B Polynomial parts 



Divergent parts can be put here 



B.l Complex or QCD 

The coefficients of polynomials part for at NNLO. 



71 = ?>2Kl^ + ?,2Kl^n - 9QK[^ - '^QKl^n + 96/s:^5 + ^2KlQn + UK^ - 64K;^ 

29LV 19L2 L2 80LL[,n 64LLI ^^^^^ 
+mKi + 32K'„n +_ + _ + _ ^ + ^ - 64LL; 

3n 3 3 n 

224LL 



r-\2 



-32LL^ - 64LL^n + ^ + 256L^L^n + 256L^L^ - 512L^L^n - 512(L^) 

2 /nV^ 164571^ _ 35 An^ _ 181 \ 
^""^^ ^ 1J28~ ~ 2^ ^ "9" ~ 

/229Ln2 4L 26L SOL^n 256L5 32L^ 368L^ 
1^^16~ ^ ~ ~9 9~ ^ 3 9~ 

_ 8L|n ^ 256L| ^ 256L^ ^ 64L|n _ ^ 192L^ 



9n 9 9 n * n / 

72 = -32X[3 - 32Kl^n + UKl^ + QAKl^n - IQK[^ - SKl^n - IQKl^ - 2,2Kl^ 

-96^3^ - 16X3^3 + 32^3^ - Y^L'n^ " ^ " ^ + + n 

ITGLL; 248LL5 ALU.n IQLU. 32LLI 20LL:n 

+ ^ + ^ + + + + ^ + 48LL^ 

3 3 3 n 3 3 

+8LLln - 32LlLln - 32{Llf + MLlLln + 64L^L^ 

/ 445Ln2 37L 16LJ; 80L: 512L^ 16L^n 

■'" +48LJ + 8LJn| 




^ 79 10368 n2 36 

73 = 2i^[i + 8i^[3 + SK[^n - 16K{^ - UKl^n + IQK^s + 48i^^' - 4i^3"i + 8K^ 

^ ^ 288 16 ° 3n 3 ^ 

8LL^n 20LL^ 



+vri6 ^ — H ^ - 31L5 ^ 

'^V 216 16 3 9n 9 ^18 

^ 56L5 20L^ 2Lgn\ 



9n 3 3 y 

n^TT^ 44571^ 37r2 23 



72 5184 16 192y 

74 = 2ir[, + 4K> + 4i^3^ + 8^5^ + 2ir,^ + 2i^> + - — - 



^ - 2LL\n H H — 

6n 6 3n 3 3 



[29Ln^ , 17 L ISU^n ^ 56L^ lOL^ 31L^ UL^n 



75 = 



76 = 




^^^^ 1^ 3456 165888 96 ^1152) 
The coefficients of polynomials part for Cp at NNLO. 



TT 



-12SKI, + 128KI, + UK; - 11^ + ^ - 1^ - 64LL> + 



•n? n 'n? n 



32L^ 192LLI 96LLI 96LL^ . . 64LL^ 



+ ^ + 64LL> ^ + ^ + 16LLI - ULLln + 



64:LL^j 192LLI 

n Ti^ 

'r\2„ -I noA T'r rr 



+ ^^^—^ - ^""^ " - 32LLI - 2^Q{Llfn - 256L^L^ + 1024L^L^n + 512L^L^ 



+512^^5^;: - 1024(Lg)'n - io24l;:l^ 

/ 12L 4L 64LS 64/.^ 64L^ 64L^ 64L^ 64L^ 192L^ 
+7ri6 T- H T' ^ \ T' + 



+^16 I — + 

-UKl^n + 128X[8 - 192/^2" + 32^20 + 64/^21 + 64is:^2^ - 32/^^2 + 64^^ 
37L2n 3L2 416LLS 96LL^ 

— T — — h ;5 — — i-b L/L/Q — 

36 n 3n^ n 



10 3 



-32^3^3 - 64^9^ + — + ^— ^ - 16LLS - 16LL> 



64LL!; 416LL^ 128LL^ 32LL!, 



T—^ + , , ' ^ + 16LL> + i + 24LL^ - 32LL^n 

6n 6n^ 3 n 

-48LL^ + 128(L^)2n + 128L^L^ - 256L^L^n - 256L^L^ 

/ 31Ln L 608LS ^^^^ 64L^ 16L^ 6O8L5 176L5 

+^^n"^"n + -9;^"^^^^^-^"-9;r+-9;^^ 9- 



32L^ 
n 



+32L> + i + 24L^ - 32L^n - 48L^ 



36 



H 7-^ + 32LL[n H ^ H H — H — ^ 



3 n 3 3n Sn^ 

AQLl 



-IQLLln - QLLl 



(llLn 224LS 8L[; ^^^^ 16L^ 8L^n 
-9;^ + -9- l^L^n 6L5 



+7r 



2 /625n 25n7r2" 
I 1296 ~ 432 



16 




5 r 2„ 

^4 = 4X[5 + 4ii'[6n + 8/^2^9 + - 4LLo - 4i^^2^ - 2^^5 

+7ri6 + 2LI - 211) + 

= -16X2^ + 2XI + 2X,^ + ^-H^-8L^n-^-iI^ 
^ ^ ^ ^ 192 3 ^33 

/lOlLn 29L[, 20L^n 97L^\ 2 ASnyr^ 115n\ 
1^92 9~ " ^^^^ " ~9 18"J + i W " 6912 J 

192" 



= 6KI - 2K; + —— - 3LLI - LLl 



Ln _ 5Ll\ 2 /^5n7r2 _ 437n\ 

— - 3Lo - — j + 7r,6 ^ - 



B.2 Real or adjoint 

The coefficients of polynomials part for Bp at NNLO. 

71 = 32KI2 + 64:Kl^n - 96K{j - 192K{sn + QGfT^g + 64X^6"' + 64i^3 - 64/r^7 

29L2n2 19L2 83L2n ITL^ igL^ 80LLSn 



+96i.;, + 64i.> + _ + _ + _- _ + _ ^ 

32LLS 8OLU0 n.rrr '^'^^LL^ arrr 32LLS 5611^ 64LL^ 



3n 3 '3 3n 3 3 

+12^ _ 48LL| ^ 64LL| _ _ _ ^ 112LL| _ 

3 n 3 n 



r\2 



+512LlLln + 2hQLlLl - 1024L^L^n - 512(L: 

f229Ln^ L 623Ln L 155L SO^n 128LS 224LS 

+^16 (^^^ + + - 7- + ^ - + 



216 n2 216 n 108 9 9n 9 

32L^ 368L^ SL^n ^ 128L^ 56L^ ^ 256L^ ^ 64L^n 32L^ 




y n 

■^^^ - 12811 - ^'^Lln + 



2^2 



^^^^ 27 ^ 1728 8n2 ^ 5184 ^ 8n ^ 54 1296 ) 
-32K[^ - 6AK[^n + 64i^[7 + 128i^[8^ - 16i^[9 - IGK^^^n - IGK^^ - 32K;^ 

- IQKl, + 32^3^, - ^IlW - - ^ + ^ - ^ + 2ALLln 
^ 36 8n2 72 8n 16 ° 

, 8LLS _^ 80LLS ^ 176LL^ ^ 248LL^ ^ ALLln ^ SLL^ _^ 32LL5 ^ 32LL^ 



1645n2 35 10763n 27 ISvr^ 2149\ 



3 3 3 3 n 3 

^20LL^ + + 48LL^ + SLLln + ULL^ - 64L^L^n - ?>2{Llf 

+128L^L^n + 64L^L^ 

, 445Ln2 317Ln 743L SLl 272LI 80L^ 512L^ 




853 \ 

72"'" 10368 ' 4n2 216 3456 4n 432 " 1296 ) 

2Kl^ + 8Kl^ + 16/s:[4n - 16/^(7 - 24i^[8n + 16K;^ + 48ir^ - AK^^ + 8/^5" 

29LV 23L2n lOlL^ ^^^^ lOLL^ 

+2K + AK^n H 8LL^n H ^ - 16LL^ 

^ ^ 288 72 96 ° 3n ° 

^ 40LL^ ^ + ^ - 6LL5 - 8LLI - LLln - 2LLI 



3 ^ 3 3n 

/49Ln2 503Ln 2867L 7L[;n 28L[; 43L:, 62L; ^ 

+7ri6 ^ \ — H - 31L; 

'^l 216 864 1728 3 9n 3 9 ^ 



25L5n 28L 



+7r 




72 5184 864 288 1728 20736^ 



2is:[i + 8i^[8n + AKl^ + 8/^5" + 2Klj + 4is:> + 



288 72 96 3 



IQLUr, ALUr, ALLl , IQLUn ALLl SLL\ LL'^.n 2LLI, 

H - 2LL\n H H — H 

3n 3 3 ^ 3n 3 3 3 3 

/29Ln2 23Ln 137L UU^n 28LI UL^ lOL^ 31L^ 
\^T6 SeT ^ 1728 9~ ^ ~9n 9 3 9~ 

llL^n 28L^ 



9n 




1 , , 421^2 riTT^ 289n IItt^ 1091 \ 



108 3456 288 10368 1728 20736 



- 8K^ - AK^ - - - — + 5LLS 5LL- 

^ ^ ^ 1152 2304 32 12 2 2 

H ^ H — H ^ 

24 8 

19Ln2 13Ln 307L 2Lln TU^ 7L[ 35L^ 5L^n 7L5~ 




+ -7^ + -r^ + — + -7^ + 



2304 768 1152 9 3 3 6 18 12 



. /nV^ 1015n2 29n7r2 10313n IQtt^ 71 



V3456 165888 3456 165888 1152 13824/ 

^ ^ 384 768 96 12 6 6 

25LL^ ^ 23LL5n hLLl 

' ^ 24 
65Ln 307L 
6912 
1933n2 

'3456 ~ 165888 ' 3456 ' 165888 
The coefficients of polynomials part for Cp at NNLO. 

5i = 128X[on - 128X[8 + 128X^ - UK^^ - 128K;^ - 2bQKl^n + 128X^6 + 384K^7n 

-1287^3^, + 128ii'Io + 64ir,^ - ^ + ^ + f + ^ - ^ - 64LL> + ^ 

4n-^ 2n 2 n"' n 

-MLLl + - l^LLl - + 64LLln - + mLU, + 

77, n.^ n. nr 





5, = 



n n n 

+8LL; - ULLln + - 32LL^ + - - 16LL^ - 512(L^)2n 



II II ir 

-256LILI + 2048LlLln + 512L^L'8 + 512L^L;; - 2048(L^)2n - 1024L^L^ 

/or T T ^RTr QOTr 1«rr QOTr 



n'' n 




16Li^32L;_ 



n 



V4n3 2n2 An 6/ 
-128X[on + 128ir[8 - 192X^ + 32X^0 + 641^^1 + l28Kl^n - 32K^^ + 6AK^^ 

-32^33 - 64X, + ^ - — - — + - 8LL, - — — + 48LL, 

-16LL> _ ^ ^ ^ ^ 1^ _ ^ ^ le^^., ^ 1^ _ ,,^^r 

6n 3 3n^ 3 n 

+12LLI - 32LLln - 2ALLI + 256{Llfn + 128L\LI - bl2L\Lln - 256L^L^ 

/ 31Ln L 5L 152Uq 32L; 81^ 8L^ 

+71-16 -; ttt + ^ " - 32L n + 32L: + — ^ 

24 4n 18 9n2 ^ n ^ 9n 9 

^^^-^ + 32L^n + - 16L^ + 12L^; - 32L^n - 24L^ j 





TT 



■2 373n 25 bn'^ 




12n 2592 16n 216 

13L2n 



^3 = 32X[on + 4X[5 + 8Xi> - 32X[8 + 96X^ - 8X2^9 + 16X^2 + 16i^9 



/llLn 115L 56L; 4LJ „„„ 8L; SL^ij 20L; 



20L5 56LS^83L|_^^^.^_3^; 



9 9n2 9 



2 f 25n7r2 625n ITtt^ 1451 \ 
^^^^ \ 86F ^ 2592 ^ "86? ^ 5184 J 



= 4X[5 + 8Kl,n + 8/^2^9 + __ + _- 2LLS - 4LL^n - LL^ 

16 48 ^ ^ 



nTT 



2 '7^ _2 



2 / iiii 7n 71^ 11 \ 
^^^'^ V 288 ^ 64 ^ 288 ^ 192/ 

^ ^ ^ 384 48 6 ^ 3 6 

/lOlLn 13L 29LJ; 20L^n 97L^ 



+7r 



V 384 384 18 ' 9 36 

2 / I9mx^ 115n TT^ 349 \ 



16 



1152 13824 ' 1152 13824 J 



6X1 - 2KI + + — - 



384 48 2 



/Ln 13L 3LS 5L5\ 
VT28 + 384 ~ ^ ~ l2"J 

2 /5n7r2 437n tt^ 349 



+7r 



16 



V1152 13824 1152 13824/ 

B.3 Pseudo-real or two-colour 

The coefficients of polynomials part for at NNLO. 

71 = 32ir[3 + 64X[4n - 96K[7 - 192ir[8n + 96X^5 + UK^qU + 64ir^ - 64X^7 

29L2n2 83L2n 191.2 ITL^ 19^^ SOLL^n 

32LL5 80LL;_ 224LL|_^^^ 3^ 

3n 3 ^ 3 ^ 3n 3 3 



+80LL^ + 512L^L^n + 256L^L^ - 1024L^L^n - 512(L^)2 

/229Ln2 623Ln L L 155L 80L[;n 128L;: 224LJ; 
+7^16 —777 ^77- + ^ + - + ^777^ ^ + + 



216 216 ^2 n 10 



32L^ 368L^ 8L5n 128L^ 
"~3 9 3~ ^ 9n 

.^_l^_128L^_32L^n + 




75 = 



+7r 



2 / 1645n2 10763n 35 27 IStt^ 2149\ 



16 



1728 5184 8n2 8n ^ 54 1296 J 



27 

■2 2 



17, o o SSL^n 3L2 SL^ 43L2 



-96K - 16KL + 32^7 - —L\^ + —=z tt^ - ^ ttt- + 2ALLln 

^ 36 72 8n2 8n 16 ° 

8LL5 80LL(; 176LL^ 248LL^ ALLln 8LLI 32LLI 32LLI 

3 3 3 n 3 3 

2QLL^n 14[ 

+ + iSLLl + SLLln - ULL^ - 64L^L^n 

/ 445in.2 317Ln 743L 8K 272L; 80L; 

" ^32- + -loT - W + - IT - V + V 

2(^22 3865n2 35n7r2 1837n 3 3 9171^ 



72 10368 216 3456 An? An 432 




73 = 2K{^ + + 16ii'[4n - 16/^1^7 - 24X[8n + 16/s:2'8 + 48^3 - AKl^ + 8/^5 

29L2n2 23L2n lOlL^ lOLL^ ^^^^^ 

+2K + 4irrn H 8LL> H ^ + 16LL^ 

^ ^ 288 72 96 ° 3n ° 

- 40LL^ - + + 6LL^ - SLLl - LL> + 2LL^5 



3 ^ 3 3n 

/49Ln2 503Ln 2867L IL^n 2SLI A^Ll. 62L\ 

+7ri6 \ — H H - 311/^ 

216 864 1728 3 9n 3 9 

25L^n 28L^ 20L^ 2^n hU\ 

— H + 3L^ ^ 2_ ^ 5 

18 9n ^ 3 3 3 / 



2 /nV^ 445n2 67n7r2 59n 18571^ 2705 \ 
^^^^ V^72~ ^ 5184 ~ 864 ^ 288 ^ 1728 ^ 20736 j 

74 = 2X[, + 8X> + 4i^3i + 8i^5^ + 2X7^ + 4X>+Hg^ + ^'" 
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C Scattering lengths 
C.l Complex or QCD case 
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C.2 Real or adjoint case 



1 _ 1 n 

1113 113 



TTCn = +X2 h 

° I 32 32n 8 



+xl 014- — as - —a2- T^TT- ai + a4 + -a3 + -a2 + —ai 

\ n 2n 8n szn 2 8 32 

111 Tl Tl Tl Tl 

0/7 1 I n 
+"^^"4- 32 + 32;?-16;^+4 + T 



3/2 1 1 1 3 ^^,1^1, 

\ n n 2n 8n 32n 2 8 

+ 32 + 2 75 + 74 + 2 73 + g 72 + 32 ^1 + 2 ^ '^s + 2 8 

71 Tl Tl^ Tl^ \ 

+ — ^1 + 2 n 74 + - 71 + 4 55 + 53 + ^ ^2 + — 5i 1 

3/7 15 3 y, 1 3^ 3^ 

''''''''' [32^^~32^^^^^'~^^'^^^'^^^' 

29 6 2 2 6 6 2 2 ^ 

_Zv _Zvg Zvg -\- Zvg + Zv^ _Zvg L2 



32n n n n n n ' n n 

6 21t7 
+ -L - 13 L^ + 3L^ -8L^- 13 L^- 13 + — L 

+16nL^ + 6nL^ -4nL^ -26nL^ + 16nL^ + 6nL^ + 6nL^ 

+16 n L + 16 + 12 + 12 LI + 16 

8 

+12 LI + 12n^ L\ + 16n^ LI - — L + 16 L^ + 16 
+16n3L^ + 16n=^L^ ) 



+''^^''^[ 48 8n3 + 4n2 + 4n 8 + 48 + 12 J 

2 2 s / 89 1 1 5 n 49n2 Sn^X ^ 



85 3 3 5 3n 2971^ 



n 



3^ 



1 1 l l n n 

— (Xa H — Pa A Pi dA 

3 24 3 24 ^ 3 24 



, /2 , 1 1 2 1 1 2n n n ' 
(^-5e + -54 + ^5,--76 - 374 -^72-y76 - 3 74 -^72^ 

^''^^"'^ (28L^^~3^^+3i^^"9i^^+2^^'3+2^ 
^ 1 2 1 2 2 

288n 3n » 9n + 9n + 9n 27n 27n ^ 

9n ^ 27n ° 1296 9 ® 9 ^ 54 " 27 ^ 
5 2 2 10 125n 4n 7n 

4n n 8n 4n 2n 7n2 

27 4 27 3 27 ^ 27 ^ 27 ° 864 27 ^ 

+ 432 

„ , / 13 1 7 281n 151n2 2bn^\ 



^ V 96 96n3 288n2 1728 2592 5184 / 
2 2 3 M73 1 _ ^ _ 1 169n _ \ 

+7r TTie^s (^10368 ^ 432n3 ~ 864^2 ~ 576n ^ 10368 ^ 432 ~ 1296 ) 



TTUq = +X2 I — - — h 



(105) 

1 1 n\ 
3"2 ~ Ten ^ 16 J 

2 113 113 



+xl a4 as - — Q!2 - —r- ai + a4 + - as + - a2 + — ai 

\ n n An Ion 2 8 32 

n ^ n n \ 

+P4 + — Pi + nai + —ai] 
16 16 J 

of 7 1 1 n 

3/^4 2 1 1 3 r 1 f n 

+2^2 75 74 73 - ^ 72 - 7^ 7l + <^4 + — 5i + 2 75 

\ n n n An Ion 16 

113 n \ 

+74 + 2 73 + g 72 + ^ 71 + ^74 + 71 I 

3 / 7 ^ 19 ^ 12 4 12 12 13 



8n3 16n2 n^ n^ n^ 16n 

12 8 4 8 12 8 8 12 

-^8 -^6 "I -^5 "I -^4 -^3 -^2 -^1 -^0 

nnnnnnnn 
+^L-13L^ + 4L^ + 3L^-4L^-13L5 + 4L^ + 4L';-13LS 

+ — L + 8 n + 8 n L;: - 2 n - 8 n + 8 n + 8 n 
8 

+SnL\ + %nLl-—L + An^Ll + An^Ll + An^Ll-—L\ 

L 41: / 



48 4n3 8n2 8n 48 16 24 
41 5 13 1 29n 
288 ^ 24^ ~ 48n2 ^ 16n ~ 288 ^ 96 ^ 



Tu) 



TTQ. 



FS 



1 



+a;2 \a3 + -a2 + —ai + ^4 + ^ Pi] + ■^'^w x] 



16 



16 



8 



+xl + ^4 + — ^1 + 4 75 + 73 + - 72 + — 7i 



16 



16 



+7ri6 X2 



L - —L + - L - ALl-SL' + iL' + SLl-iL' 



^8n? An ' 8 
-8L^-8L^ -4L5 



+7^16 



7ra 



MA 



17 1 1 n 
" 24 ~ 4^ ^ 2^ ~ 24 
1 1 / 



13 



TT 



+ 



1 



24 



+ 



144 24n2 
1 



12n 



7n 
144 



288 288n2 



1 



1 



72n 



9n 



9n 



2 
27ri 



9;?^^+27;?^^+ 



27n2 
1 



— + L--LI 
9n ^ 2592 9 ^ 



2 2 1 

9 ^ 9 ^ 27 ^ 



27 



9^9° 



7n 
1296 



17 



1 1 n 

^^■^2 I 2592 96n3 ~ 144n2 ^ 648, 



Wnlexl - 



1 



+ 



1 



1 



1 



n 



1296 432n3 864n2 864n 5184, 



7ra, 



MS 



+4 



1 



1 



1 



1 



+X2 U4 + — 5i - 2 75 - - 73 - - 72 - — 7i 



16 



8 



32 



+7ri6 X2 



16n 32 ^ 



An 



ALl-Ll 



+ALI + ALI-LI 



3 Pseudo-real or two-colour case 



1 1 n 
111 



vrao = +a:. --- — + - 



+X9 I aA 0:3 ao ai — ad as ao Oii 

^\ n 2n 8n 32n 2 8 32 

111 Tl Tl Tl Tl 

+A - jft - g/?2 + ^ A + 2llQ4 + - Q, - - A - -ft - ^ft 

2/7 1 1 n 
+vri6X,(^-- + ^ + — -- + -J 

3/^2 1 1 1 3 

\ n n 2n 8n szn 2 8 

1 1 3 ^ r n ^ n . 
+ ^'^i-2 75 -74 --73 --72 - — 7i-2n55--(53--52 

-32 + 2^74 + g 71 + ^5 + + — (52 + — (5i j 

3 / 7 ^ 15 ^ 3 1 ^ 3 3 29 ^ 

+^-^H^'^^32;^^^^'^^"^'^^^^'^^^^'^°"32;^'^ 

6 27-r 2rri^rri^rr 2;^ 2rri^rr 

+ ~ 1^8 1^6 -^5 "I -^4 "I -^3 -^2 -^1 "I -^0 

nnnnnnnn 

55 21n 
-— L-13L^ + 3L^ + 8L^-13L^-13LS+— L-16nL: 
32 lb 

+6n LI + An LI - 26n LI - 16n LI + 6n LI + 6n LI - 16n LI 
L + 16n2 - 12 LI - 12 + 16n^ L^ - 12 L^ 

-12n'^L{ + 16n'^Uo- —L + 16n^Ll + 16n^Ll + 16n^Ll 

+16 L\^ 

2 3 /85 _ 3 _ 3 5 _ 3n _ 29n2 n^N 
+^16 2^2 [^--^-^ + ^-Y~^^12) 

2 2 3 89 1 1 5 n 49^2 5n=^\ , , 

2/211 3 113 

+^4-n"^-n"^-4;^"^-16;^"^-"^-2"^-8"^-32"^ 



r 



1 n \ 

16 16 J 

3/^4 2 1 1 3 r 1 f ^ 

+2^2 75 74 73 - 72 - 7^ 7i + (^4 + — ()i - 2 75 - 74 

\ n n n An Ibn 16 

113 n \ 

-r^-8^^-32^^ + ^^^+16^V 

3 / 7 ^ 19 ^ 12 4 12 12 13 ^ 

12 8 4 8 12 8 S r'- _L ■'■^ 

i Ivg ijg i H Ivg -^1 -^0 

nnnnnnnn 
L - 13 - 4 + 3 + 4 - 13 - 4 - 4 - 13 

+ — L-8nL^ + 8nL^ + 2nL^-8nL^-8nL^ + 8nL^ 

+8nL^ -8nLj; + y L + 4n2L^ + 4n2L^ + 4n2L[; ^ 
'37 5 13 3 31n Sn"^ 7n^\ 



2 2 3/41 5 13 _ 1 _ 29n _ Sn^X 
''^^''^V 288 + 24;^ + 48;?"T6;^" 288 "% + 144j' ^^^^^ 

/l 1 1 1 n n\ 

1,3"^+ 24 "^+3^^+ 24 ^^"3"^" 24 "V 
2/11 1 1 7n 

1, ~ 864 + 288;? + 288;; + 864 ~ 432 J 

0/2^112 1 1 2n n 
+X2 (^3^6 + 3^4 + ^52+376+374 + ^72-^76-374 

n \ 
-24^^j 



5 ^ 1 ^ », 2 ^„ 1 ^„ 2 ^„ 5 ^„ 2 



r 



L + — L";,- — L"^- — LI + — LI + — LI - — L2 



288n 3n ^ 9n ° 9n ° 9n ^ 27n 27n 

9n ^ 27n ° 1296 9 ® 9 ^ 54 ^ 27 ^ 
5 2 2 10 12577 477 7n 

-^^3 + 9^2-^^i-^^o + ^^ Y^6 + ^^5 

4n 77 877 477 277 7772 772 

"27^^ 27^^ 27^^ 27^^+27^° 864 ^ 27^^ 



+7r^g X2 



432 
I 96 



281n 151^2 



rr 



96n3 288n2 
173 1 



+ 



+ 



1728 
5 



2592 
1 



25^3 \ 
5184; 
169n 



+ 



n 



10368 432n3 864^2 576n 10368 432 



1296 ) ' 



7ra, 



FA 



1 

H Xo 

16 

- q;3 - - q;2 - — ai + P4 + — Pi + -TTie 



1 



1 



+^2 U4 + — ^1 - 4 75 - 73 - - 72 - — 7i 



16 



16 



8 



877,2 4^ 
+8L^ + 8L^ -4L[; 

1 



n 



17 ^ 

24 ^ 4^ ^ 2^ ~ 24 ^ ^ 



7n 



13 



7ra 



MA 



+xl 



+xl 



1 



144 144 

-±-\ 

288 ■ 288^2 i 



24n2 12n 



3^6 + 354 + ^52 



+7ri6 x^ 



V 288^3 
1 ^ 



9n 



72n2 

i:;; + 



1115 

+ T. — ?r -^s — T. — ?r -^^ + — ~ Lr, 



2>'n? 



9n2 



27n2 



9n 



27n 



9n 



11 
2592 



27n2 
1 , 



72n 

2 2 1 2 2 1 7n 

T r I T r r r r r r r T ^ _L r 



9 



9 



+7rlexl 



9 
17 
2592 

1 



+ 



+ 



27 
1 

96ri3 
1 



+ 



27 
1 



+ 



9 



n 



+ 



144^2 648 J 
1 1 



nan — 



2 ^1296 432^3 864^2 
1 



864n 



9 



5184 i 



1296 



32 



X2 



1 



1 



1 



+X2 2 "3 + g «2 + ^ «1 + ^4 + A + ^TTie 



+xU + ^4 + — 5i + 2 75 + - 73 + - 72 + — 71 



16 



32 



V 

J_ _ J_ _ J_ _ lln 2/j_ J. 1_ 7n\] 

_Sn 96 8n2 96^^ V 576 ^ 48^2 48n^ 576yJ' ^ ^ 

D Loop integrals 
D.l One-loop integrals 

We use dimensional regularization here throughout in d dimensions with d — A — 2e. 
We need integrals with one, two and three propagators in principle. These one 
propagator integral is 

A{m^)^- [ -^—1—. (116) 

The two propagator integrals we encountered are 

d'^q 1 
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B{mi,m2,p ) ^ - JTT 

I J (ZTT 



{27rY{q^-ml)iiq-py-ml)' 

1 f d'^q q,, 



{27iY {q^ - ml){{q - p)^ - ml) 
= p,B,{mlmlp'), 
r> ( 2 2 N _ 1 [ _d^ q^ 

= Pu.Pi^B2i{mj,ml,p'^) + g^^B22{mj,ml,p'^) . (117) 

All the cases with three propagator integrals that show up can be absorbed into 
the two-propagator ones by moving to the real masses rather than the lowest order 
masses. This provides a consistency check on the calculations. 

The explicit expressions are well known 



2 ^2 „2N _ 1 mllnjiiif) - inlHiiil) \ o 2 2. , x 

B[m^,m2,p ) - Y^l^o ml -ml j + 'J {^i^^2,P ) + (-'[^) , 

167r^ Jo \ mfx + m^^l — x) J 

C = ln(47r) + 1 — 7 and Aq = 1/e + C. The function J{m\^m\^p'^) develops an 
imaginary part for p^ > (mi + 7712)^. Using A = mf — m\, E = m\ + m| and 
1/2 — p4 _|_ ^4 _|_ ^4 _ 2p'^m\ — 2p^m\ — 2m\m\ it is given by 

9v -/ 9 9 9x /A S\ m? z/ , (r)2 + z/)2 - 

(32.^)jK,™i,/) = 2+(-- + -lln-|--l„l^j-i-^.(119) 



The two-propagator integrals can all be reduced to B and A via 
Bi{ml,ml,p^) = -^[A{ml)-A{rril) + {ml-ml-p^)B{ml,ml,p^)), 

B22{mlmlp^) = ^^^^-—^[A{ml) + 2mlB{mlml,p^) 

+ (m2 -m^ -p )Bi{m^,m2,p )) , 
S2iK,m^,p2) = ^(A{ml) + mlB{ml,ml,p^)-dB22{ml,ml,p^)) . (120) 

The basic method used here is the one from Passarino and Veltman 
D.2 Sunset integrals 

Sunset integrals have been treated in many places already, in general and for various 
special cases. We use here a method that is a hybrid of various other approaches. 
We only cite the hterature actually used. We define 
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H{mlmlml;p^) = ((1)), 

Hf,^{ml,ml,ml;p'^) = 

= p^,p„H2i{mj,ml,ml;p'^) + g,j,„H22{ml,ml,ml;p'^) . 

(122) 

By redefining momenta the others can be simply related to the above three: 

{{r^)) = Pi,Hi{ml,ml,ml;p'^) , 
{{r^ir„)) = p^,p^H2i{ml,ml,ml;p'^) + g,^„H22{ml,ml,ml;p'^) , 

{{%r^)) = P^iPvH2z{m\,ml,ml]p^)+ g^^H2A{rn\,ml,ml]p^) , (123) 

with 
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+H22{mlmlml,p^). (124) 

The first two follow from interchanging q and r and the third from the fact that it 
is proportional to g^i, or p^Pi,, which are both symmetric in /i and u. The last one is 
deribed using 

(g^r^ + ri^q^) = (g^ + - p^,) (g^ + - p^) - q^q^ - r^,r^ - p^,p^ 

+2p^ {q, + r,) + 2p, {q^ + r^) (125) 



and redefining momenta and masses on the r.h.s.. In addition we have the relation 

2tt / 2 2 2 2\ I jTj / 2 2 2 2\ 

p H2i[mi,m2,mo/,p ) + dH22{'m^,m2,m^;p ) = 

mlH{ml, ml, ml; p") + A{ml)A{ml) . (126) 

which allows to express H22 in a simple way in terms of H2\ . There is also a relation 
between Hi and H 
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Hi{mi,m2,m-i;p ) + iii(m2, m^, mg;^ ) + Hi[m^,mi,m^\p ) = 

H{mlmlmlp'') , (127) 

which allows to write i7i(m^, m^,m^]p^) = 1/3 H{m^, m^, m^; p^) in the case of equal 
masses. The function H is fully symmetric in mj, m| and m\, while ifi, if2i and H22 
are symmetric under the interchange of m\ and mg. 

We only need the sunset integrals at p^ = m\ = m\ = m\ and their derivatives 
w.r.t. p . These have been calculated using the methods of [29]. With Hid = -^H^ 
we obtain 
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Hi and i/i,^ follow immmediately using (127). 



D.3 Vertex integrals 

The vertex diagram (16) in Fig. 2 is the most difficult two- loop diagram in (fxj) 
scattering, and it can also appear in other process. The two loop integral for the 
equal mass case can be written as 

i"^ J J (r^ — m^) • [(r — q)^ — m^] ■ (s^ — m^) • [(s + r — p)^ - m^] 

(134) 



The Lorentz decompositions of the vertex integrals are [34] 

((1)) = V, 
{{fn)) = PnVii + q^yi2, 

((Sa*)) = Pm^13 + ?m^14, 

{{rtiSy)) = gf,uV25 + ^26 + qu.quV27 + qnPuV28 + P^tquV29 , 

( (S^S,.)) = 9^^V210 + + g;ig;.14l2 + {qnPi. + Pnqv)V21Z , 

{{ri^r^ra)) = {g^^vPa + Q^iaPv + gvaP^i)Vu + {g^ivqa + g^iaqv + gvaqn)Vz2 

+Pu.PuPaV33 + qi^quqaVzA 

+{Pf.Puqa + P/i^i^Pa + q^iPvPa)Vz^ + (g^g^^Pa + g^^Pi^ga + Pt,quqa)Vm , 

{{r,j.rvSa)) = gfMuPaVsj + gfMuqaVss + (gf^aPu + guaP^,)V3Q + {g^^aqy + S'l/a?^) V'siO 

+Pt,PvPcyi,ll + g^gj^gaVais + Pi^PuqaV^l^ + g/.gi.PaV'aM 

+ (PA'gi^ + qi^Pu)PaV3l5 + (p^gi. + g^Pi.)gaV3i6 , 
((^^/xSySa)) = P/.5'j^aV3i7 + q^^guaVsiS + (fl'/.i^Pa + 5'a»qP!/) ^319 + (fi'/.^.^a + S'/xagi/) ^320 
+P/iPi.PaV321 + g^gi^ga 1^322 + P^iquqcy32-i + gMP'^Pa^324 
+PlJ.{Piyqa + qvPa)V32b + g/.b^'^a + ?:^Pa)^326 • 

(135) 

The {{s^SySg)) does not show up in 00 scattering. Most of those Vi functions have 
been calculated analytically in [20] except the {{s^Sy)) and {{r^SySa))- We have 
calculated the rest of them in this work, which are V210 — V213 and V317 — V326- Again, 
the methods of [29] were used here, somewhat extended to the cases at hand. We 
have compared our results with the numerical evaluation for general masses described 
in [34]. The quantity B'^ is the next term in the expansion of B in (118) but these 
terms always cancel in the final result. 
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Where the J and function are defined as 
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The J and ki vanish at s = and are well behaved for s — )■ cxo. They have disconti- 
nuities in the derivative at threshold but there no poles there. The functions ki are 
constructed using the arguments and methods of [35]. 

All the ki above show up at intermediate stages of the calulations but in the 
final result k^{s), . . . , kg{s) always appear multiplied by powers of s and can thus be 
removed. 

Finally, to get the scattering lengths we need to expand these functions around 
t,M = and s = 4. The expansion using s = 4(1 + g^) around s = 4 reads up to 
order q^. 
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